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Lyapunov’s Linearization Method 


Define A = ae) ‘ and decompose f(x) as 


x= 


Ig(x)| 


f(x) = Ax+ g(x) where "a 


— Oas |x| > 0. 


Theorem: The origin is asymptotically stable if R{A;(A)} < 0 for 
each eigenvalue, and unstable if R{A;(A)} > 0 for some eigenvalue. 


Note: We can conclude only local asymptotic stability from this lin- 
earization. Inconclusive if A has eigenvalues on the imaginary axis. 


Proof: Find P = P’ > 0 such that ATP + PA = —Q <0. Use V(x) = 


x'Px as a Lyapunov function for the nonlinear system x = Ax +(x). 


V(x) = x1 P(Ax + 9(x)) + (Ax + g(x))"Px 
=x"(PA+ A'P)x + 2x! Pg(x) 
< —x7Qx + 2|x|||Piilg(x)| 


Amin(Q) |x? < xT Qx < Amax(Q) |x? 


V(x) S —Amin(Q)|x/? + 2I|Plilxllg(x)| 


Is(x)| 
|x| 


|x] <r => |g(x)| < |x|; see the illustration below for the case x € R. 


Since 


— Oas x — 0, for any y > 0 we can find r > 0 such that 


¥|x| 
> |x| 


Thus, |x| < r(y) = V(x) < —Amin(Q)|x/? + 27||PIllx/?. 


Choose 7 < Sue so that V is negative definite in a ball of radius 


r(y) around the origin, and appeal to Lyapunov’s Stability Theorem 
(Lecture 8) to conclude (local) asymptotic stability. 
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Region of Attraction 


Ry ={x: p(t,x) > 0} (1) 
“Quantifies” local asymptotic stability. Global asymptotic stability: 
Ry = R". 


Proposition: If x = 0 is asymptotically stable, then its region of at- 
traction is an open, connected, invariant set. Moreover, the boundary 
is formed by trajectories. 


Example: van der Pol system in reverse time: 


x1 = —X2 

. 3 (2) 
Xq = Xy — XQ + X5 

The boundary is the (unstable) limit cycle. Trajectories starting within 

the limit cycle converge to the origin. 


Example: bistable switch: 


Xy = -ax, + X92 

(3) 
; xy 3 
Xo = — bx 
2 14x 2 


x2 


> X1 


Estimating the Region of Attraction with a Lyapunov Function 


Suppose V(x) < 0 in D — {0}. The level sets of V inside D are 
invariant and trajectories starting in them converge to the origin. 
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Therefore we can use the largest levet set of V that fits into D as an 
(under)approximation of the region of attraction. 


A 


D tei Vx cy oC Ry 


This estimate depends on the choice of Lyapunov function. A simple 
(but often conservative) choice is: V(x) = x’ Px where P is selected 
for the linearization (see p.1). 


Time-Varying Systems 
Khalil (Sec. 4.5), Sastry (Sec. 5.2) 
t=f(tx) f(t,0) =0 (4) 
To simplify the definitions of stability and asymptotic stability for 
the equilibrium x = 0, we first define a class of functions known as 
"comparison functions." 


Comparison Functions 


Definition: A continuous function a : [0,00) — [0,00) is class-K if 
it is zero at zero and strictly increasing. It is class-K’. if, in addition, 
a(r) + coasr > om, 


A continuous function B : [0,00) x [0,00) — [0,00) is class-KL if: 
1. B(-,s) is class-K for every fixed s, 


2. B(r,-) is decreasing and B(r,s) — 0 as s — oo, for every fixed r. 


Example: a(r) = tan~!(r) is class-K, a(r) = r°,c > O is class-Ko, 
B(r,s) =re-® is class-KL. 


Proposition: If V(-) is positive definite, then we can find class-K 
functions «1(-) and #2(-) such that 


ai (|x|) < V(x) < a2(|x|). (5) 
If V(-) is radially unbounded, we can choose @;(-) to be class-Ko. 
Example: V(x) =x'Px P=P™>0 


ay(|x|) = Amin(P) |x|? a2(|x|) = Amax(P)|xI?. 
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Stability Definitions 


Definition: x = 0 is stable if for every € > 0 and to, there exists 6 > 0 
such that 
|x(to)| < d(to,€e) => |x(t)| <e Vt> fo. 


If the same 6 works for all to, i.e. 6 = d(e), then x = 0 is uniformly stable. 


It is easier to define uniform stability and uniform asymptotic stabil- 
ity using comparison functions: 


¢ x = (is uniformly stable if there exists a class-K function a(-) and 
a constant c > 0 such that 


|x(#)| < a(|x(to)|) 


for all t > to and for every initial condition such that |x(to)| <c. 


¢ uniformly asymptotically stable if there exists a class-K’L B(-,-) s.t. 


[x(t)| < B(|x(to)|,t — to) 
for all t > to and for every initial condition such that |x(to)| <c. 


¢ globally uniformly asymptotically stable if c = oo. 


* uniformly exponentially stable if B(r,s) = kre~*S for some k, A > 0: 


|x(E)| < klx(to) |e AO) 
for all t > to and for every initial condition such that |x(to)| <c. 


Example: Consider the following system, defined for tf > —1: 


= =e 
1+t 


x(t) = x(tp)elo >” = x(ty)eloB(+s)l° 

1+ to 

1+t 

|x(t)| < |x(to)| == the origin is uniformly stable with a(r) = r. 


log Ltto 
= X(ip)e°S HF = x(fp) 


The origin is also asymptotically stable, but not uniformly, because 
the convergence rate depends on to: 


1+ to x(to 
x(t) = x(to) — = (Fo) 


increasing to 


t—to 


